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Abstract. Even though the first momenta i.e. the ensemble average quantities in canonical 
ensemble (CE) give the grand canonical (GC) results in large multiplicity limit, the fluctuations 
involving second momenta do not respect this asymptotic behaviour. Instead, the asymptotics 
are strikingly different, giving a new handle in study of statistical particle number fluctuations 
in relativistic nuclear reactions. Here we study the analytical large volume asymptotics to 
general case of multispccics hadron gas carrying fixed baryon number, strangeness and electric 
charge. By means of Monte Carlo simulations we have also studied the general multiplicity 
probability distributions taking into account the decay chains of resonance states. 

1. Introduction 

The surprising success of the statistical hadronization model in describing the 
multihadron production in high-energy nucleus-nucleus reactions (see e.g. [1] and 
references therein) as well as in elementary particle collisions [2] gives rise to a question 
if the fluctuations, especially the second moments of multiplicity distributions, follow 
the statistical resonance gas results as well as the mean multiplicities (the first 
momenta of multiplicities) . Even though the CE first momenta reach the GC values 
quite fast with increasing multiplicities [3] , we have found that the variance of these 
values always stay different in the CE [4, 5] due to the lack of charge bath, i.e. due to 
the lack of fluctuations of the relevant net charges. 

In section 2 we pose and solve the theoretical problems involved, and give results 
for the general hadron gas. The section 3 is devoted to the numerical studies employing 
the methods developed in previous section, and to the comparison of those with the 
ones given by the Monte Carlo (MC) procedure explained in Ref. [6]. 

2. Theoretical Studies 

The central quantity in our discussion is the scaled variance of the multiplicity N: 

This is a finite measurable parameter even at the infinite volume limit, and we will 
show it losing its volume dependence in large systems. It is worth stressing that for the 
Poissonian distribution corresponding to Boltzmann GC ensemble, u> = 1. Although 
we keep the theoretical discussion on the level of Boltzmann approximation for brevity, 
we show also the full quantum statistics results in the section 3. 
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In the CE, the partition function does not factorize to one-species parts, but 
the different species are bound by the group factors giving delta functions over 
corresponding exact conservation laws. Recall the form of the partition function for 
canonical hadron gas carrying baryon number B, strangeness S and electric charge Q 

[7]: 
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The first and second momenta of a given multiplicity of a hadron subset h can be 
calculated by inserting fictitious fugacity to the generating function Zgc(T, V, {Aq s — > 
e «0»}) m equation (2) and performing the derivatives 
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Using these, the scaled variance can be written as the sum of Poissonian part and 
the canonical correction: 



UJ h = l + 
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For the large systems, the canonical correction factors can be solved by means of 
the asymptotic expansion, where the saddlcpoint of the f(x) (3) coincides with the 
values of GC fugacities [8]. In the general case of multi- hadron gas carrying total 
charges {Qk} = {B, S, Q) the asymptotic expansion of the partition function involves 
diagonalization of the Hessian matrix 
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Because this is real and symmetric, the diagonalization is done by the orthogonal 
matrix A : H' = diag(/ii, hi, /13) = A T HA, where hk are the Hessian eigenvalues, 
and A is constructed of corresponding eigenvectors. Applying these and the notation 
{N%)gc f° r GC average multiplicities, the large volume limit of the Eq. (5) is found 
to be 
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Figure 1. LEFT: Scaled variances in the canonical (c.e.) neutral pion gas as functions of GC 
one-particle partition functions (i.e. GC multiplicities). GC Poissonian value is denoted by 
g.c.e. 

RIGHT: Scaled variances in the thermodynamical limit in nuclcon-pion gas (T = 120 MeV) 
as functions of baryochcmical potential fis- uqb stands for QB canonical result, whereas for 
u)q the GC approximation is used for baryon number conservation. In the former case, /ib is 
taken only as a saddle point in the projection integral (3), not as a Lagrangian coefficient. 



where the summations of Qk,i are over relevant charges B, S and Q, and the A l Q k are 
the elements of the transformation matrix A T . 

As an easily evaluable approximation, we consider the scaled variance for the 
number of positively (h + ) and negatively (h~) charged particles as well as the 
sum of all charged particles neglecting the multi-charged particles and the mixing 
between different charges (non-diagonals in the Hessian). The relevant Hessian term 

is Hqq = X Q 2 Q + J2j z i^jQj(Qj - 1) j which leads to H QQ = X Q 2 [Q + 2(/i~)gc] = 

A Q 2 K^ + )gc + (h-) G c}- This yields the terms Q i Q j /((h+) GC + (/i")gc) inside the 
Qk, m summations in the Eq. (7). Applying this to equation (7) yields 
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v-oc * (h+) GC + (h-)ac v^oo \{h+) GC + {h-)Gc) K ' 

This result serves as a "rule of thumb" for any charge under consideration - it applies 
trivially to baryons and strange particles too - but in quantitative comparisons one 
must relax the approximations used above. In case of one-charge system, such as pion 
gas, Eq. (8) is an exact result for thermodynamical limit. One sees immediately, that 
for a neutral system limy^oo ui± = 0.5 instead of GC limit u± = 1. When charge 
density grows indefinitely, the limiting values are limy^oo uj + = 0, limy^oo w_ = 1 
and limy^oo ui c h = 0. 
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3. Numerical Results 

In some special cases, Eq. (5) can be easily evaluated for any volume [4, 5], not 
only in the asymptotic limit. Thus, in Fig. 1 we quote the neutral pion gas charged 
particle fluctuations as functions of GC multiplicities [4]. Also in Fig. 1, we show the 
asymptotic large volume results as functions of baryochemical potential in the nucleon- 
pion gas in i3Q-canonical calculation [5] together with the approximation where the 
B is handled in GC manner, see Eq. (8). 
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Figure 2. Full quantum hadron gas results in T = 160 MeV, Q/B = 0.4. LEFT: Scaled 
variances for charged, positive and negative hadrons as functions of baryon density. The lines 
depict the asymtotic results, and dots are for MC results together with statistical errorbars. 
The arrows show the change of oj_ and u> c h, when the strong decay chains are taken into 
account in addition to primary production. 

RIGHT: Scaled variances for primary baryons, antibaryons, strange hadrons and for antistrange 
hadrons. 

In Fig. 2 we turn to the full -BSQ-canonical ui calculations as functions of baryon 
density for temperature T = 160 MeV. The isospin condition is chosen to correspond to 
the PbPb and AuAu reactions: Q/B = 0.4. We take into account all the resonances up 
to the mass ~ 1.9 GeV. The left panel shows the asymtotic results for scaled variances 
of primary negative, positive and charged hadrons together with Monte Carlo results 
for V = 200 fm 3 , which is checked to be large enough to give asymptotic results. Also 
shown are the values for oj_ and w c h after the strong decay chain. In the right panel 
we present the scaled variances for primary baryons, antibaryons, strange hadrons and 
anti-strange hadrons. 
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